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Abstract. The effect of poloidal asymmetry of impurities on impurity transport 
driven by electrostatic turbulence in tokamak plasmas is analyzed. It is found that in 
the presence of in-out asymmetric impurity populations the zero-flux impurity density 
gradient (the so-called peaking factor) is significantly reduced. A sign change in 
the impurity flux may occur if the asymmetry is sufficiently large. This may be a 
contributing reason for the observed outward convection of impurities in the presence 
of radio frequency heating. The effect of in-out asymmetry is most significant in regions 
with low temperature gradients. In the trapped electron mode dominated case also an 
up-down asymmetry can lead to a sign change in the peaking factor from positive to 
negative. The effect of ion parallel compressibility on the peaking factor is significant, 
and leads to positive peaking factors in regions with high temperature gradients, even 
in the presence of in-out asymmetry. 
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1. Introduction 

Accumulation of impurities in the core of fusion devices would have detrimental effect 
on fusion reactivity due to increased radiation losses and plasma dilution. Significant 
effort has therefore been spent during recent years to identify plasma conditions in which 
accumulation can be avoided. One of the most promising methods for obtaining flat or 
hollow impurity density profiles is to use radio frequency heating. This has been shown 
to work well in various experiments [H El El lU El El tl] but the physical mechanism by 
which the change of the direction of the impurity convective velocity occurs has not yet 
been clearly identified in spite of the various efforts that have been made [El El ED]- In 
particular the reason for the flat impurity density profiles in ion cyclotron resonance 
heated (ICRH) discharges in JET has been debated for many years [U El E] • In a recent 
paper [TT] it was shown that a poloidal asymmetry could lead to a significant reduction 
of the impurity zero-flux density gradient (also called the peaking factor), and even a 
sign change in the impurity flux, if the asymmetry is sufficiently large. This suggests a 
possible explanation for the avoidance of accumulation of high-Z impurities with ICRH. 

In this paper a numerical investigation of the effect of poloidal asymmetry on 
impurity transport is presented. In particular, the dependence of the impurity density 
peaking factor on charge number, inverse ion- and electron temperature scale lengths 
and inverse electron density scale length is analyzed. Also, the role of parallel ion 
compressibility is studied. The results are benchmarked to GYRO [12] in the poloidally 
symmetric case. 

It is found that inboard accumulation gives rise to negative peaking factors (outward 
impurity convection and hollow impurity profile), both in ion temperature gradient 
(ITG) and trapped electron (TE) mode driven turbulence. The sign and magnitude of 
the peaking factor will be shown to be sensitive not only to the asymmetry strength 
but also to the temperature gradient, and it is more pronounced in the low temperature 
gradient region characteristic for the plasma core. The peaking factor is expected to be 
negative for moderate to high impurity charge numbers (above Z ~ 15) for JET-like 
experimental parameters [1]. In TE mode dominated cases up-down asymmetry can also 
give rise to negative peaking factors if the electron temperature gradient is sufficiently 
large. As noted in previous work [8], in the poloidally symmetric case, that parallel 
compressibility has a significant effect on the peaking factor. As we will see in this 
paper this effect is even more pronounced for asymmetric impurity densities, and will 
lead to negative peaking factors for lower asymmetry strengths than without taking into 
account parallel compressibility. 

The remainder of the paper is organized as follows. In Sec. [2] we describe the 
mechanism behind the impurity poloidal asymmetry that arises in the presence of ICRH. 
In Sec. Elthe model for calculating the quasilinear impurity flux and the peaking factor in 
the presence of poloidal asymmetry is presented. In Sec. @] the parametric dependences 
of the peaking factor are analyzed by presenting scans over relevant parameters such as 
charge number, and temperature and density scale lengths. Also, the importance of the 
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impurity parallel compressibility is demonstrated. Finally, the results are discussed and 
summarized in Sec. [5j 

2. Poloidal asymmetry 

Poloidal impurity asymmetries in tokamaks can arise due to various reasons: e.g. 
difference in impurity source location, toroidal rotation or neoclassical effects. There 
is a wealth of experimental evidence for poloidal asymmetries [131 EH HSl EH EH EE] • 
In the plasma core in-out asymmetries can arise due to the presence of radio frequency 
(RF) heating (henceforth "in-out" and "out-in" asymmetries will refer to the situations 
when the maximum of the poloidally varying impurity density is at the inboard and 
outboard sides of the plasma, respectively). A detailed physical explanation of why 
ICRH favors inboard accumulation, together with the description of the experimental 
results, is given in Ref. [13]. The underlying principle is that the asymmetry is a result 
of the increase of the hydrogen-minority density on the outboard side. These particles 
tend to be trapped on the outside of the torus and the turning points of their orbits 
drift towards the resonance layer due to the heating. The poloidal asymmetry in the 
hydrogen-minority density gives rise to an electric field that pushes the other ion species 
to the inboard side. In the case of highly-charged impurities, this effect is amplified by 
their higher charge Z. 

The RF-induced accumulation of minority ions on the outboard side leads to a 
corresponding impurity accumulation on the inboard side by the following mechanism. 
If the plasma consists of electrons, bulk ions, impurity ions and RF-heated minority 
ions then it can be expected that all species except the minority ions are Boltzmann 
distributed (the dynamics of the minority ions is strongly affected by the heating). If 
a particle species a is Boltzmann distributed, the poloidal variation of the density is 
na/ n a,o — —G a 4>/T a , where the tilde denotes the variation on the flux surface, e a is the 
charge and T a is the temperature of the species. Quasineutrality requires that 




n e0 1^1 + YJ ~ u do ^1 - jr J - Zn zo ^1 - -j^ 



where the subscript zero indicates the density where the potential <fi vanishes. If <fi is 
normalized so that nno + n m + Zn zo — n e0 = 0, the poloidal variation of the impurity 
density becomes 

n z Ze<f) _ Zn H /n D0 

n Z0 ~ T z ~ 1 + (Tj/T e ) + (n Z0 Z 2 /n D0 ) ' 1 j 

Since the poloidal variation in the impurity density has the opposite sign to that of 
the minority ions, the accumulation of the latter on the outboard side gives rise to an 
electric field that pushes the other ion species to the inboard side. Simulations of the 
hydrogen ion distribution function in the presence of RF heating with the Monte Carlo 
code FIDO [19] show that a considerable out-in asymmetry in the hydrogen ion density 
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can be expected, which is sufficient to account for the observed in-out asymmetry in the 
impurity density 

In the tokamak edge, where the plasma is sufficiently collisional, also steep radial 
pressure or temperature gradients can give rise to an in-out asymmetry. These effects 
have been observed in e.g. Alcator C-Mod [H], and it has been shown that the 
observations are in qualitative agreement with neoclassical theory [201 12H 1221 [23] . The 
sign and magnitude of these asymmetries depend sensitively and nonlinearly on magnetic 
geometry, fraction of impurities in the plasma and rotation. Neoclassical theory also 
predicts an up-down asymmetry, which is caused by the ion-impurity friction. 



3. Impurity flux 

Since impurity transport is usually dominated by drift-wave turbulence, in this work 
we focus on the effect of the impurity poloidal asymmetry on impurity transport driven 
by microinstabilities. We assume that the processes that cause the asymmetry are 
not affected significantly by the fact that the cross-field transport is dominated by 
fluctuations. For simplicity we consider only the collisionless, electrostatic case. The 
quasilinear impurity particle flux is given by 

r , = -|lmM*], (2) 

where Im[-] denotes imaginary part, kg is the poloidal wave-number, n z is the perturbed 
impurity density, 0* is the complex conjugate of perturbed electrostatic potential <j). 

The perturbed impurity density response in an axisymmetric, large aspect ratio 
torus with circular magnetic surfaces can be obtained from the linearized gyrokinetic 
equation [21]. In the absence of collisions and in the limit v\\/qR{u — wdz) 1, the 
gyrokinetic equation can be solved iteratively to find the non-adiabatic part of the 
perturbed impurity distribution 

g =(l- Vl1 9 9 \ ZC ^ z0 - - U * z JJz U 

\ qR(u - lo Dz ) 86 qR(u - uj Dz ) 86 J T z to - lo Dz 

where 6 is the extended poloidal angle, f a0 — n a / (y/irvTa) 3 exp(— x 2 a ) is the equilibrium 
Maxwellian distribution function, x a = v/vT a is the velocity normalized to the 
thermal speed vx a = (2T a /m a ) 1//2 , n a and m a are the density and mass of species a, 
uJ* a = -k e T a /e a BL na is the diamagnetic frequency, uj a = w* a 1 + [x 2 a - §J L na /L Ta , 
L na = — [<9(lnn a )/(9r] _1 , L^ a = — [<9(lnT a )/9r] _1 , are the density and temperature 
scale lengths, una = — ke (y]_/2 + v jj) D(8)/u ca R is the magnetic drift frequency, 
D{6) = (cos 9 + s^sin^), u ca = e a B/m a is the cyclotron frequency, B is the equilibrium 
magnetic field, q is the safety factor, s = (r/q)(dq/dr) is the magnetic shear, r and R are 
the minor and major radii, Jo is the Bessel function of the first kind and z a = k±v±/u ca . 
Including the Boltzmann part of the distribution, the perturbed ion density response 
becomes [M 125] 

v\\ 8 v\\ 8 



n z Zecj) r 3 Zef z0 J (z z ) 



, d 6 v- 

n z T z J n z T z 
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(3) 

where parallel compressibility is represented by the term containing the two 9- 
derivatives. 

The zero-flux impurity density gradient (the peaking factor) can be obtained 
from (T z ) = with the perturbed impurity density taken from (|3]). Here (• • •} = 
(l/27r) fl n {- ■ -)d0. If the impurity density varies over the flux surface, in general we can 
write n z = n z oV(9), where V(9) = 1 if the impurities are evenly distributed over the 
flux surface. A poloidal asymmetry can be modelled by V(9) = cos 2 (J2) > wnere ^ 
is the angular position where the impurity density has its maximum. This asymmetry 
function was used also in previous work [TTJ[T5]. Figure [U shows the asymmetry function 
V(9) as a function of normalized 9/tt, for n = 4 and different values of S. 




Figure 1. V{ff) as a function of normalized 9/tt, for different values of the peaking 
angle 6. The asymmetry strength is n — 4. Solid line (red) represents out-in, dashed 
line (blue) is up-down and dash-dotted line (black) is in-out asymmetry. 



The impurity flux can be calculated numerically, by solving the velocity-space 
integrals in the perturbed impurity density ([3]) in the expression for the impurity flux 
(j2j), without the constant energy resonance approximation or the assumption on the 
smallness of the finite Larmor radius parameter. Then the impurity peaking factor 
ajL? nz can be obtained by setting the impurity flux to zero, or (lm[n z 0*]) = 0. Here, 
the fluctuating density is given by Eq. (jHJ), a is the outermost minor radius, and a/L° nz 
is calculated as 

a/Ll = p&m, (4) 



where 
S n (9) 
and 
S T (9) 



J d 3 x z e x2 *Jq[z; 
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with the operator representing parallel compressibility 

x 4 ® x z\\ 9 

9 = R(9)(oo - u Dz (9)) d9R{9){uj - oo Dz (9)) d§' ( ) 

where R{9) = Rq(1 + ecos9) with e = r/Rg, 5 P = 2a 2 mi/ (q 2 m z T z ), r z = T e /T z , 



b = (p s /Z)y2Az /r z Ai, p s is the ion sound Larmor radius, u* z = —kep s /Zr z , udz = 
-2k ePs (a/R)(x 2 z± /2 + x 2 4 )D(9)/Zt z , z z (9) = x z± bk ± (9), and k ± = k e Vl + s 2 9 2 . All 
frequencies marked with bars are in c s /a units, where c s is the ion sound speed, and the 
main ion and impurity temperature gradients are assumed to be equal a/L Tz = a/L Ti . 



4. Parametric dependences of the peaking factor 

In the calculations presented in this section we have used the following local profile 
and magnetic geometry parameters: r/a = 0.3, R/a = 3, kgp s = 0.3, q = 1.7, 
a/L ne = 1.3, a/L m = 0.86, T t /T e = 0.85, a/L Te = 0.2, a/L Ti = 0.3, s = 0.22 
and p s /a = 0.0035. These parameters have been taken from a typical JET-shot 
with ICRH [1]. The impurities are assumed to be present in trace quantities, in 
the sense that Zn z /n e <ti 1 (n z /n e = 10~ 3 is used in the simulations). This is the 
baseline case in our study, and these parameters will be used unless otherwise stated. 
The perturbed electrostatic potential and eigenvalues have been obtained by linear 
electrostatic gyrokinetic calculations with gyro . 



4-1. Charge number dependence 

The peaking factor as a function of charge number for various impurity poloidal 
asymmetries is shown in Fig. In the poloidally symmetric case (solid line in Fig. Ek), 
the peaking factor is not sensitive to the charge number, as has been noted before, in 
both fluid and gyrokinetic simulations of ITG turbulence dominated transport, without 
taking into account the poloidal impurity asymmetries [H [26] . The situation is similar 
also in case of up-down asymmetric impurity populations, while the peaking factor 
increases slightly in case of outboard accumulation. However, in agreement with the 
conclusion of [11], the black dotted line in Fig. |2^i shows that impurities experience 
outward convection (corresponding to negative peaking factor) if the impurity density is 
accumulated on the inboard side. The change in the peaking factor becomes stronger as 
the asymmetry is increased, as it is illustrated in Fig. |2]d. Note, that the strength of the 
asymmetry is also expected to depend on Z, and usually it is larger for heavy impurities, 
as it was shown in Sec. [2j According to Eq. (prj), in the limit of trace impurities, the 
poloidal variation of the impurity density is proportional to the charge number. 

In order to examine the importance of parallel compressibility in determining the 
charge dependence of the impurity peaking factor we separate the terms independent 
and proportional to 5 P (representing the parallel compressibility) in Eq. fl5]) as S n (9) = 
Sk(9) + S*>(e), where 

SM = V{9)m f d 3 x z e-^J 2 [z z (9)] ^ (8) 
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Figure 2. Peaking factor as a function of charge number for different values of peaking 
angle 5 (a) and asymmetry strength n (b). The rest of the parameters are taken from 
the baseline case, for which the eigenvalue is tu = (0.001 + i 0.068)c s /a. In both 
figures the solid line represents the case of poloidally symmetric impurity distribution, 
(a): n = 4 - out- in asymmetry (dashed), up-down asymmetry (dash-dotted), in-out 
asymmetry (dotted), (b): in-out asymmetry - n = 2 (dashed), n = 4 (dash-dotted). 



sew = s, [*x**JMQ\M, ~ v VM>-y* m , 0) 

J UJ — uj Dz {V) 

and in Eq. © as S T (9) = S^(9) + S P T \6), where 

Sh(e) = V(6 )m [ d* Xz e^jZ[z z (e) f " °*- {xl : V ( % a/LTZ , (10) 

J UJ—UDz\V) 

Using the above notations the peaking factor can be rewritten as: 

n/r o (MgWI) + (MSFW]) n?] 

1 nz (Jm[Si(9)<t>*]) + (Im[S£°(0)^]) ' 1 J 

Figure [3] shows these four expressions as functions of impurity charge for the 
symmetric and in-out asymmetric cases with an asymmetry strength of n = 4. In 
the denominator, the terms independent of S p , i.e. S^(6), are the dominant contributors 
in both the symmetric and asymmetric cases. These are proportional to 1/Z. But in the 
numerator, the balance between the terms proportional and independent of 5 P is very 
different between the symmetric and asymmetric cases; in the symmetric case S^(6) 
is the dominant term, while in the asymmetric case S^(6) is dominant. This latter 
term, which is due to parallel compressibility, changes sign at moderate charge numbers 
(Z ~ 15), which in turn leads to a sign-change in the peaking factor. These results show 
that if an in-out asymmetry is present, the parallel compressibility effects become more 
important than the other terms, and therefore, have to be taken into account. 

In the following we will concentrate on the peaking factor for nickel which was 
one of the impurities studied in Ref . [1] . Figure H] shows the peaking factor for nickel 
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z z 



Figure 3. Terms proportional (dashed lines) and independent of S p (solid lines) 
in Eq. (|12[) as functions of impurity charge Z. (a) and (c) show the terms in the 
denominator and (b) and (d) show the terms in the numerator, (a) and (b) correspond 
to the poloidally symmetric case, and (c) and (d) correspond to in-out asymmetry with 
strength n = 4. 



for various accumulation maxima 5 and asymmetry strengths n. The peaking factor 
increases slightly for up-down accumulation but the more dramatic change - including 
a sign-change - is expected only for inboard accumulation. The sign change occurs when 
n ~ 2. In Ref. pXJ it was shown that without taking into account the effect of parallel 
compressibility the sign change occurs when n ~ 3. This is in agreement with the results 
obtained here when parallel compressibility effects are neglected, i.e. by setting 5 P = 
in Eqs. ([5]) and ([6]) (as will be shown later in Fig. [T3|) . 

Assuming that by increasing the RF heating and therefore the temperature 
gradients, the asymmetry will increase as well, e.g. n oc a/L^ the results shown in 
Fig. can be linked to the experimental results illustrated in Figs. 6 and 7 in Ref. [I]. 
In this reference the experimental observation of the effect of ICRH heating on the nickel 
transport in the plasma core at JET is discussed. It is shown that the application of 
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Figure 4. Peaking factor for nickel as a function of n (a) and 8 (b). (a): 
symmetric impurity density (solid, red), up-down asymmetry (dashed, blue), and in- 
out asymmetry (dash-dotted, black), (b): symmetric impurity density (solid, red), 
n = 2 (dashed, blue), n = 4 (dash-dotted, black). 



the ICRH leads to an outward impurity flux (negative peaking factor). As the applied 
heating is increased and therefore, the temperature gradient peaked, the outward flux 
is further increased leading to more negative peaking factors. 

4-2. Temperature gradient dependence 

The eigenvalues and electrostatic potentials as functions of ion- and electron temperature 
gradients are shown in Figs. EHH1 As expected, if we increase the ion temperature 




Figure 5. Real and imaginary parts of u> = uj r + i"f as function of a/Lxt (a) and 
a/Lxe (b) obtained by GYRO for the baseline case. Blue lines (with circle markers) 
represent the real part, blue lines (triangle markers) correspond to the imaginary part 
of the eigenvalue. The frequencies are normalized to c s /a. 



gradient, the turbulence becomes more ITG-dominated (the real part of the mode 
frequency u r is negative), while if we increase the electron temperature gradient, TE- 
mode driven turbulence will dominate (u r is positive). The shape of the imaginary 
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part of the potential lm[0] varies strongly by increasing the temperature gradient for 
both the ITG and TE-mode dominated cases, while the real parts of the potential Re[0] 
are not modified significantly, see Fig. EJ The imaginary part of the potential plays 




Figure 6. Real and imaginary parts of the electrostatic potentials in the cases 
corresponding to Fig.[5]for different ion- (a,b) and electron (c,d) temperature gradients, 
(a) and (c) show the real parts and (b) and (d) show the imaginary parts of cf> 
(note the different scale). Red lines (circle markers) show the temperature gradient 
corresponding to the baseline case and blue lines (triangle markers) show the case with 
the larger temperature gradient. 

an important role in the parallel compressibility terms and therefore the change in the 
temperature gradient will modify the impurity peaking factor considerably. 

Previous works highlighted the difference in peaking factors in ITG and TE 
dominated cases, and concluded that ITG-dominated turbulence will always generate 
inward pinch of impurities, while in TE-mode driven turbulence outward convection 
in the plasma core (for r/a ~ 0.2) can be expected. Both linear [7J |8] and non- 
linear [27] gyrokinetic simulations have shown that the latter is due to the contribution 
from the parallel dynamics which can reverse the impurity convection from inwards 
to outwards for modes propagating in the electron direction. These results are in 
agreement, under some conditions, with the experimental observation that the impurity 
convection changes sign from inward to outward when a strong central peaking of the 
electron temperature arises as a response to strong localized central electron heating. 
However, these results cannot explain the outward convection of impurities observed in 
experiments with RF heating where the ITG-mode is the dominant instability. 

Our results show, that the peaking factor is negative in both ITG- and TE- 
dominated cases if the impurities accumulate on the inboard side. However, the peaking 
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factor is influenced very strongly by increasing temperature gradients and will become 
positive again above a certain temperature gradient. 




Figure 7. Peaking factor for nickel as a function of ion temperature gradient for 
different values of S (a) and n (b). In both figures the solid line represents the case 
of poloidally symmetric impurity distribution; in figure (a) this case is compared to 
GYRO simulations (black diamonds), (a): n = 4 - out-in asymmetry (dashed), up- 
down asymmetry (dash-dotted), in-out asymmetry (dotted), (b): in-out asymmetry - 
n = 2 (dashed), n = 4 (dash-dotted). 



Figure [7] shows the ion temperature-gradient scaling of the peaking factor and the 
diamond symbols show these values obtained by GYRO simulations in the poloidally 
symmetric limit, which present good agreement with our results. The peaking factor 
is not very sensitive to the ion temperature gradient as long as the impurity density 
is poloidally symmetric or if it is up-down or out-in asymmetric. However, if the 
impurity density is in-out asymmetric then the sign of the peaking factor is negative 
for low temperature gradients and positive (and even larger than in the poloidally 
symmetric case) for large temperature gradients. The threshold for the sign change 
for the experimental scenario studied in this paper is a/ L^i — 1. As the strength of the 
asymmetry, n, is increased for an in-out asymmetry the modifications of the peaking 
factor becomes stronger, see Fig. [7b. 

As the electron temperature gradient is increased the turbulence becomes more TE 
dominated. Figure M shows that the absolute value of the peaking factor is reduced up 
to a certain temperature gradient and then starts to increase again, regardless of the 
sign of the asymmetry. The inboard accumulation gives negative peaking factors also in 
this case. One of the most interesting differences compared to the ITG-dominated case 
is that also an up-down asymmetry can lead to a sign change for a/LiTe — 0.5. This 
can be a contributing reason for the observed flat impurity density profiles in plasmas 
with electron cyclotron resonance heating (ECRH), where the in-out accumulation 
mechanism described in [13] and in Sec. [2] would not be present. Up-down asymmetries 
have been observed in electron cyclotron (EC) heated plasmas [15]. The physical 
mechanism for these up-down asymmetries was attributed to neoclassical effects, but 
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since the magnitude of the asymmetries reported in Ref. [TH] was very large it is possible 
that there are other reasons that may be related to the ECRH as well. 




Figure 8. Peaking factor for nickel as a function of electron temperature gradient 
for different values of 8 (a) and n (b). In both figures the solid line represents the case 
of poloidally symmetric impurity distribution; in figure (a) this case is compared to 
GYRO simulations (black diamonds), (a): n = 4 - out-in asymmetry (dashed, blue), 
up-down asymmetry (dash-dotted, green), in-out asymmetry (dotted, black), (b): in- 
out asymmetry - n = 2 (dashed, blue), n = 4 (dash-dotted, black). 



4-3. Density gradient dependence 

The density gradient scaling of the peaking factor is shown in Fig. [9j The peaking factor 
is slightly increasing in the case of poloidal symmetry and in the cases of up-down or 
out-in asymmetries, as the electron density peaking is increased. Also here, the in-out 
asymmetric impurity density leads to negative peaking factor if the strength of the 
asymmetry is sufficient. The peaking factor is fairly insensitive to the density gradient 
in this case. It is interesting to note, that even though the ITG-mode is the dominant 
instability here, (see Fig. [TO]) , the peaking factor remains negative as the ITG-mode 
is becoming stronger in the case of an in-out asymmetry. This can be contrasted to 
the case when the stronger ion-temperature gradient led to a sign change and positive 
peaking factor, see Fig. [7k. The difference between the two cases can be understood by 
comparing the shape of the imaginary part of the electrostatic potential lm[0] in Fig. [TTb 
with that in Fig. [6b- It can be noted that by changing the temperature gradient, lm[0] 
changes significantly, while changing the density profile does not result in a significant 
difference. Note that it is mainly the part of the potential which is close to 6 = tt which 
is important, and that is considerably different for the two temperature gradients shown 
in Fig. [6b and therefore also the result for the peaking factor changes dramatically. The 
difference in the electrostatic potentials is the underlying reason for the difference in the 
trends seen in Figs. [7J and [9] as the ITG-mode become stronger by increasing the a/Lxi 
or a/L ne , respectively. 



Possible mechanism responsible for observed impurity outward flow under radio frequency heatingYi 




Figure 9. Peaking factor for nickel as a function of electron density gradient for 
different values of S (a) and n (b). In both figures the solid line represents the case 
of poloidally symmetric impurity distribution; in figure (a) this case is compared to 
GYRO simulations (black diamonds), (a): n = 4 — out-in asymmetry (dashed, blue), 
up-down asymmetry (dash-dotted, green), in-out asymmetry (dotted, black), (b): in- 
out asymmetry - n = 2 (dashed, blue), n = 4 (dash-dotted, black). 




Figure 10. Real and imaginary parts of the eigenvalues as function of a/L ne obtained 
by GYRO . Red lines (with circle markers) represent the real part, blue lines (triangle 
markers) correspond to the imaginary part of the eigenvalue. The frequencies are 
normalized to c s /a. 

4-4- Effect of parallel compressibility 

It has been shown previously that when the transport is TE-mode dominated parallel 
compressibility effects generate an outward contribution to impurity anomalous flux 
which can, under certain plasma conditions, cancel out the inward contributions, leading 
to zero or even negative impurity peaking factor [TJ |S]. In this subsection, we examine 
this effect by neglecting the parallel compressibility terms, i.e. terms proportional to S p 
in Eqs. (jSJ) and (Q. Note that in this case only the absolute value of the potential enters 
in the expressions, while in the case with parallel compressibility both the imaginary 
and real parts and their derivatives are important. Figure [T2l illustrates the square of 
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Figure 11. Real- (a) and imaginary (b) parts of the electrostatic potentials for two 
different density gradients. 



the electrostatic potential |0| 2 for different values of a/Lnpi and a/L ne . We note that 




Figure 12. The square of the absolute value of the electrostatic potentials, \4>\ 2 for 
different ion temperature gradients and density gradients. 

the absolute value of the potential is not significantly different in the various cases, 
so therefore we do not expect dramatic dependence on the temperature and density 
gradients. Figure [T3l shows the peaking factor for various S and n. The effect of poloidal 
asymmetry in this limit is in agreement with our previous work in Ref. [H] and also, 
is similar to that in the previous sections where the parallel compressibility effects are 
considered. As seen in Fig. [T3k . in the absence of parallel compressibility effects an 
in-out asymmetry will lead to a negative peaking factor (outward impurity flux). An 
increase of the poloidal in-out asymmetry will increase the outward flux of impurities 
as shown in Fig. [T3b . Note that in the case with parallel compressibility, the sign 
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change in the peaking factor occurs for broader range of S and for lower asymmetry 
strength (compare Figs. 0] and [13]). Also an up-down asymmetry can lead to a slight 
reduction of the impurity peaking factor. This is different from the case where the 
parallel compressibility was taken into account, see Fig. [4k. 
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Figure 13. Peaking factor for nickel as a function of n (a) and S (b), without taking 
into account parallel compressibility, (a): symmetric impurity density (solid, red), 
up-down asymmetry (dashed, blue), and in-out asymmetry (dash-dotted, black), (b): 
in-out asymmetry - symmetric impurity density (solid, red), n = 2 (dashed, blue), 
n = 4 (dash-dotted, black). 



Figure HH shows the density- and temperature-gradient and charge number scans 
for the peaking factor without parallel compressibility. The diamonds represent the 
values of peaking factor obtained by GYRO (without parallel compressibility) which 
show agreement with our results in the symmetric limit. From Fig. [14] it is clear that 
without parallel compressibility the effect of a poloidal asymmetry is almost insensitive 
to the increase of the gradients of temperatures and density or the impurity charge. 
This is mainly due to the fact that without parallel compressibility only the absolute 
value of the potential matters and that is rather similar in the various cases. 

5. Discussion and conclusions 

For tokamak operation the avoidance of central impurity accumulation is a key issue. 
Conditions in which the convective impurity flux is directed outward are particularly 
interesting. In order to find such conditions many experiments have been devoted to 
explore various techniques. One way to expel the impurities from the plasma core is to 
maintain sawtooth crashes in a controlled way by applying central ICRH [28] [29]. It 
is observed that this method will indeed remove the impurities from the plasma core. 
However, it was shown that even though sawtooth crashes hamper the accumulation 
of the impurities their contribution is less relevant compared to the effect of the ICRH 
itself (5j. Another technique which has been successful in removing the impurities from 
the plasma core and is routinely used in tokamak experiments such as in the ASDEX-U 
tokamak is the application of a very localized central ECRH [30]. Simulations with 
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Figure 14. Peaking factor for nickel as a function of ion temperature gradient (a), 
electron temperature gradient (b) electron density gradient (c) and impurity charge Z 
(d) for different values of 6, without taking into account parallel compressibility. Solid 
(red) line represents the case of poloidally symmetric impurity distribution; in figures 
(a-c) this is compared to gyro simulations (black diamonds). Out-in asymmetry 
(dashed, blue) up-down asymmetry (dash-dotted, green), in-out asymmetry (dotted, 
black) . 



both linear and non-linear gyrokinetic models have shown that under these conditions 
the electron temperature gradient is strongly peaked and therefore, modes propagating 
in the electron diamagnetic drift direction are the dominant instability responsible for 
the turbulence driven transport [27]. The interaction between the related electrostatic 
potential fluctuations and the parallel dynamics of the impurities leads to an outward 
convection of impurities. These results are in agreement with experimental observation 
in ASDEX-U for the very central region (r/a ~ 0.2) [7J. In other tokamak experiments 
for example at JET the application of the central RF heating has also been explored with 
success. In these experiments it has been observed that the dominant instabilities are not 
TE-modes but rather modes directed in the ion diamagnetic drift direction, ITG-modes, 
which theoretically should result in an inward impurity flux. It has been debated that in 
the very central region of plasmas at JET the transport is mostly driven by neoclassical 
effects and by applying a strong central heating as the temperature gradient peaks the 
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neoclassical temperature screening effects become dominant and result in an outward 
directed impurity ffux. Under some plasma conditions these effects may be the reason 
for the observed behavior [2] however, usually the observed transport is an order of 
magnitude larger than the neoclassical predictions implying that the impurity transport 
is turbulence driven [27] . 

In the present work we assume that the impurity density is in the form of 
n z(8) = n z®P{9) where the poloidal dependence is represented by V(6) = cos 2 f^f^J , 
and discuss the impact of such a dependence on the impurity peaking factor in the 
core of tokamak plasmas. Various mechanisms may give rise to a poloidal asymmetry of 
impurity density: difference in impurity source location, toroidal rotation or neoclassical 
effects. Among these is an in-out poloidal asymmetry observed in plasmas where RF 
heating is applied. The mechanism responsible for this behavior was explained through 
RF-induced accumulation of minority ions on the outboard side of the torus giving rise 
to a corresponding impurity accumulation on the inboard side as was discussed in Sec. 2. 
The strength of the impurity accumulation depends on the impurity charge (higher for 
heavier impurities) and the plasma parameters such as temperature gradient. However, 
the exact form of these dependences is not yet known and further analysis is needed in 
this area. 

We have found that an in-out poloidal asymmetry of the impurity density (5 = it) 
can lead to an outward impurity flux (negative peaking factor) in both ITG and TE 
mode dominated cases, and this effect becomes stronger as the asymmetry strength 
(n) increases. In this case the peaking factor is strongly charge dependent and it 
becomes more negative for heavier impurities. In TE mode dominated case an up-down 
asymmetry can also generate a negative peaking factor. This can be a contributing 
reason for the observed flat impurity density profiles in plasmas with ECRH, and up- 
down asymmetries have indeed been observed in EC heated plasmas [15], although the 
physical mechanism for these asymmetries is not entirely understood. 

The reason for the sign change of the impurity peaking factor in the presence of 
a poloidal asymmetry is attributed to the interaction between the poloidal variation of 
the related electrostatic potential and the poloidal dependence of the impurity density. 
If parallel compressibility effects are taken into account the imaginary part of the 
electrostatic potential is the determining factor in the sign change of the impurity 
peaking factor. 

In summary, our results suggest that poloidal asymmetries can significantly alter 
the turbulence driven impurity transport and therefore, have to be taken into account. 
These asymmetries may have a significant role in determining the impurity accumulation 
properties in plasmas with radio frequency heating. Therefore, there is a strong need 
for development of new tools in order to detect poloidal asymmetries and determine the 
effect of RF heating in their generation and the dependence of the asymmetry function 
to various plasma parameters. 
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